A lattice glass model with no tendency to crystallize 
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We study a lattice model with two body interactions that reproduces in three-dimensions many 
features of structural glasses, like cage effect and vanishing diffusivity. While having a crystalline 
state at low temperatures, it does not crystallize when quenched, even at the slowest cooling rate 
used, which makes it suitable to study the glass transition. We study the model on the Bethe lattice 
as well, and find a scenario typical of p-spin models, as in the Biroli Mezard model. 



The glass transition is a phenomenon that is observed 
when a liquid is cooled well below its melting tempera- 
ture, avoiding the crystallization of the system [|@]. The 
motion of the molecules is slowed down more and more 
by lowering the temperature, until at some point it be- 
comes so slow that the liquid appears as a disordered 
solid, with the molecules vibrating around some equilib- 
rium positions, but not diffusing anymore. Two funda- 
mental approaches have been put forward in the study 
of this phenomenon, from the point of view of dynamics 
H , and thermodynamics || . A further insight in the un- 
derstanding of the glass transition was achieved through 
the study of the p-spin glass || , a mean field model with 
p-body interactions and quenched disorder, which repro- 
duces many features of glass forming liquids. More re- 
cently, Biroli and Mezard have introduced a very inter- 
esting lattice glass model (6) , where occupation variables 
interact via p-body potentials. The mean field solution 
on the Bethe lattice shows a scenario identical to p-spin 
models, and three dimension numerical simulations show 
a behavior typical of glass forming liquids. In order to 
reduce the tendency to crystallize, the authors consid- 
ered a mixture of two types of particles. Further study 
of the model was done in j?j . Another model with p-body 
interactions and without quenched disorder, that shows 
a similar behavior, was studied in ||. 

In this paper we study a simple lattice gas model 
with only two body interactions. The model reproduces 
in mean field again the scenario typical of p-spin mod- 
els, and three dimension numerical simulations show the 
same behavior typical of glass forming liquids, including 
cage effect and two step relaxation. Moreover, while hav- 
ing a crystalline phase at low temperatures, the model 
does not crystallize when cooled, even in the pure case 
(only one type of particles) and at extremely low cooling 
rates. This makes the model more suitable to study the 
glass transition, and allows to obtain very clear results. 

The model is defined as follows. We partition the space 
in regular cells, such that not more than one particle can 
have its center of mass inside the cell. The position and 
orientation of the particle inside the cell is represented 
by a coarse grained discrete internal degree of freedom, 
that can assume a finite number q of states, and the 



interaction between nearest neighbor particles depends 
explicitly from their position inside the cell. The model 
is therefore described by the following Hamiltonian 



H = ^ n i n j < l>ij( a 'i) aj) - M y"/' 



(1) 



where n, = 0, f whether the i-th cell is occupied by a par- 
ticle or not, Oi = 1, . . . , q represents the internal position 
of the particle, and </>y (<7i, <%/) is the interaction energy 
between two particles in cells i and j with internal posi- 
tions &i and (Tj. 

It is clear that, by choosing a sufficiently large number 
q of internal positions, and an opportune interaction ma- 
trix <pij (<7i, <jj), one can approximate as closely as desired 
any model defined in the continuum, like for example a 
Lennard-Jones liquid. On the other hand, it is plausible 
that a few number of internal states may be enough to 
catch the fundamental characteristics of dynamics and 
thermodynamics of glass forming liquids. 

Here, we study a particularly simple realization of the 
model ([!]). In two dimensions, we partition the space in 
square cells, and subdivide each cell into four internal po- 
sitions. When a cell is occupied by a particle in a given 
position, a hard-core repulsion forbids the presence of 
another particle in some of the internal positions of the 
neighboring cells (see Fig. |l|). Therefore in this case 
q = 4, and the interaction (j>ij(ai, <jj) is zero if the posi- 
tions <ii and Oj are "compatible" , infinite otherwise. The 
extension to three dimensions is straightforward. In that 
case one partitions the space in cubic cells, and considers 
six internal positions instead of four. For every spatial 
dimension d, the model has a crystalline ground state 
with density 2d/(2d+ 1), on lattices with sides multiples 
of 2d + 1 lattice spacings. For example, on a cubic lat- 
tice with periodic boundary conditions it can be found as 
follows: consider the cell with coordinates (x,y,z), and 
evaluate the number a — (x + 2y + 3z mod 7): if a = 0, 
leave the cell empty; if a = 1,2,3 put a particle in the 
negative x, y or z direction, respectively; if a = 4, 5, 6 put 
a particle in the positive z, y or x direction, respectively. 

We have simulated the model in three dimensions, by 
means of grand-canonical and canonical dynamics. A 
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compression experiment on a system of size 28 3 is shown 
in Fig. ^. A simple Monte Carlo grand-canonical (vari- 
able density) dynamics is performed, and temperature 
is slowly decreased, starting from some high value. The 
results are shown as open circles in Fig. 0, for various 
cooling rates. The final high-density state is strongly 
dependent from the cooling rate, as observed in glass 
forming liquids. On the other hand, no tendency to crys- 
tallization is observed: even at the slowest cooling rate 
(T/T = — 1CP 7 ) no transition to the crystalline state was 
observed. We have also computed the density in the crys- 
talline state, starting from the ground state at very low 
temperature, and heating up slowly (diamonds in Fig. 
||). This ordered state becomes unstable at T ~ 0.150 ^i, 
where the system falls into the liquid phase. 

After having equilibrated the system at a certain den- 
sity in the liquid phase, we have switched to a canoni- 
cal (fixed density) dynamics. These simulations are per- 
formed on lattices of size 15 3 . The mean square displace- 
ment (r 2 (t)} of the particles is show in Fig. |^a. This is 
defined taking in account also the position of the particle 
inside the cell, giving to each internal position a shift 1/4 
of a lattice spacing with respect to the center of the cell. 
In Fig. [|b instead we plot the self-overlap, defined as 



Zext an d %inl b e the partition functions of the branch, 
restricted respectively to configurations where site i is 
empty, occupied on the "external" position, and occupied 
on one of the k "internal" positions (see Fig. ||b). We 
define the "local fields" ai and hi acting on the site i by 
the relations = = {Z® + Z^)/Z^. 

The branch can be seen as the result of the merging of 
the k branches ending on the neighbor sites j = 1, . . . , k, 
which leads to a recursion relations for the local fields, 
and a free energy shift: 



(«(*)> 



l]T(n 4 (tH (0)^(^.^(0)) 



(2) 



where <7i(t) are unit length vectors, pointing in one of 
the six coordinate directions, representing the position 
of the particle inside the cell, and N is the number of 
particles. As it happens in glass forming liquids, at high 
density a plateau develops both in the mean square dis- 
placement and in the self-overlap. This is the signa- 
ture of the so called "cage effect" . Particles are trapped 
in cages formed by neighboring particles, and vibrate 
rapidly around some equilibrium position, giving rise to 
the first decay in correlation functions. Only after a long 
time they manage to escape from their cages, giving rise 
to the second final decay to equilibrium. While the time 
of the first step remains finite, the time of the slow relax- 
ation diverges when the glass transition is approached. In 
Fig. ^a it is shown the diffusivity as a function of p — p Cl 
with p c = 0.844. It can be fitted by a power law, with an 
exponent 7 = 3.45. In Fig. Eb it is shown the dynamical 
susceptibility x(t) — -W((<z(*r) _ (^W) 2 ); which shows a 
behavior typical of p-spin models, and observed also in a 
Lennard Jones binary mixture J9|, with a peak at finite 
time that grows when approaching the glass transition, 
signaling the presence of dynamical heterogeneities. 

We have then studied the model on the Bethc lattice, 
namely a random lattice with fixed connectivity |||l2],|l3| . 
Each site of the lattice is subdivided in k + 1 internal 
positions, and is connected to k + 1 randomly chosen 
neighbors (see Fig. [s]a,) . For large number of sites A, the 
lattice is locally treelike, but has loops of length of the 
order log A. Consider a "branch" ending on site i, that is 

(i) 

a graph in which site i has only k neighbors, and let Zq , 
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To evaluate the total free energy, one has to consider also 
the free energy shifts when merging k + 1 branches on a 
new site, and when merging two branches, that are given 
respectively by 



k+l 



fc+1 



-f3AF { 



= l[(l + e^) + e^J2I[( 1 + e0hj ) ( 4a ) 
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-/3AF< 2 > = j + e pai + e l3a 2 + e f3(h 1 + h 2 ) 



(4b) 



The replica symmetric solution of the problem cor- 
responds to the case in which there is only one pure 
state, and the fields do not fluctuate from site to site. 
In this case, the total free energy is given by F — 
AFW - (k + l)AF( 2 )/2, from which one can find the 
specific volume and the entropy per site, that are shown 
as solid lines in Figs. || and 7] (in the case k = 5). As 
it happens in the model studied by Biroli and Mezard 
||, the entropy becomes negative for temperatures be- 
low some threshold. A similar behavior is observed for 
every k > 1. One can also find a (replica symmetric) 
crystalline solution, where the fields do not fluctuate on 
the single site, but are different on different sites, and 
the merging of the branches is done preserving the crys- 
talline structure. For k = 5, this solution appears below 
T ms = 0.1633/1, and becomes stable below the melting 
temperature T m — 0.1444/^. The corresponding volume 
and entropy are shown as dashed lines in Figs. || and |7[ 

We have then looked for a solution at the level of 1- 
step replica symmetry breaking (RSB) [ [il)||Tl|]l^ , |il| ] . In 
this case, many pure states exist, and the local fields 
fluctuate also on the single site. We work in the so called 
"factorized case" , in which the probability distribution of 
the local fields P(di, hi) on the single site is the same on 
all the sites of the lattice. The self-consistency equation, 
when merging k branches on a site i, now reads 
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P(a, h) oc J J^J dcij dhj P(a,j, hj 



xS(a- ai)6(h- hi)exp(-f3mAF) (5) 
where a%, hi and Ai* 1 are functions of Oj and /ij via the 



equations 
energy is 



), to is a real parameter, and the total free 
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One then has to maximize the free energy with respect 
to the parameter to, in the interval < to < 1. 

At temperature T = the solution of (||) is easily 
found. Indeed, in that case the fields at and hi can only 
have the finite set of values 



hi = 



(l-fc)M,-. 



0,/i 



ifai < /i 
if Qj = fi 



(7) 



The distribution P{a, /i) is then a sum oik + 2 delta func- 
tions, and the integral equation (|5|) gives rise to k + 2 al- 
gebraic equations for the weights of the delta functions. 
One has to give a finite value to the parameter ftm/i, and 
maximize the free energy with respect to it. For k = 5, 
the maximum is F — —0.8265 /i, for Bmfi — 8.078. At 
temperature T > the solution of (||) can be found it- 
eratively. We discretized the distribution P(a, h) over a 
domain of the plane (a,h), using a fine grid with spac- 
ings da = dh — /x/1024, and starting from some initial 
distribution applied iteratively equation (||) until the pro- 
cedure converged. Of course one must be careful that the 
chosen domain in the plane (a, h) covers all the support 
of P(a,h). To locate the maximum of the free energy 
as a function of m, it is useful to evaluate explicitly the 
derivative dF/dm 12 1. The results for the volume and 
entropy, and for the parameter m, are shown as open cir- 
cles in Figs. ^| and in the case k = 5. The 1-step RSB 
solution appears below Tjj ~ 0.105 fx, and becomes stable 
below Tk — 0.087 (i, where to = 1. Note that below Tk 
the configurational entropy vanishes, and therefore the 
residual entropy, at low temperatures in the 1-step RSB 
solution, has to be interpreted as a "vibrational entropy" , 
due to the "rattling" of the particles inside their cages. 
This vibrational entropy is absent in the crystalline state. 

In conclusion, we have studied a lattice glass model 
with two body interactions. On the Bethe lattice, it 
shows the same scenario of discontinuous spin glasses, 
and of other recently studied lattice models of glass. 
Moreover, it shows no tendency to crystallize when sim- 
ulated in three dimension with a simple Monte Carlo dy- 
namics, which avoids the need to introduce a mixture, as 
in many models of glass. 



After this work was done, we learned about the work 
of Weigt and Hartmann, who have studied a lattice glass 
model that shows a similar behavior |l4| . 
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FIG. 1. The model in two dimensions: the space is parti- 
tioned in square cells, and each cell can be occupied by at most 
one particle in one of four positions (little circles) . A particle 
in a given position (big shaded circle) forbids the presence of 
another particle in the positions colored in black. 
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FIG. 2. Specific volume (inverse density) as a function of 
the temperature, for a three-dimensional system of size 28 3 . 
Circles: the system is cooled starting from high temperature, 
for cooling rates (from top to bottom) T/T — — 1CP 4 , — 1CP 5 , 
— 10 -7 . Diamonds: the system is heated starting from the 
crystalline ground state, with heating rate T/T — 10~ 7 . 





FIG. 5. a) The model on the Bethe lattice (here with 
k = 3): each site is subdivided in k+1 positions (little circles), 
and connected to k + 1 randomly chosen neighbors. A particle 
in a given position (big shaded circle) forbids the presence of 
another particle in the positions colored in black, b) Merging 
of the branches j = 1, . . . , k onto the site i: the "external" po- 
sitions are those colored in black, the others are the "internal" 
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FIG. 3. a) Mean square displacement (r 2 (r)) and b) self 
overlap (q(t)) as a function of time, in a three-dimensional 
system of size 15 3 , for densities p — 0.55, 0.6, 0.65, 0.7, 0.75, 
0.76, 0.78, 0.8, 0.81, 0.815, 0.818, 0.821. 




FIG. 4. a) Diffusivity D as a function of p c —p, with crit- 
ical density p c = 0.844. The fitting function is a power law 
D — a\p — p c | 7 , with 7 = 3.45. b) Dynamical susceptibility 
x{t) as a function of time, for densities p = 0.76, 0.78, 0.8, 
0.818. 
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FIG. 6. Specific volume (inverse density) as a function of 
the temperature, on the Bethe lattice with k — 5. Solid line: 
replica symmetric liquid phase; broken line: crystalline phase; 
open circles: glassy phase in the 1-step RSB approximation. 
Arrows mark the various temperatures cited in the text. 
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FIG. 7. The same as in Fig. [|, but for the entropy per site. 
Inset: parameter m as a function of the temperature. The 
line is a guide for the eye. 
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